ANALYSIS OF THE SINGULAR SOLUTION BRANCH OF A 
PRESCRIBED MEAN CURVATURE EQUATION WITH SINGULAR 
NONLINEARITY MODELING A MEMS CAPACITOR 
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Abstract. The existence and multiplicity of solutions to a quasilinear, elliptic partial differential 
equation (PDE) with singular non-linearity is analyzed. The PDE is a recently derived variant of a 
canonical model used in the modeling of Micro-Electro Mechanical Systems (MEMS). It is observed 
that the bifurcation curve of solutions terminates at single dead-end point, beyond which no classical 
solutions exist. A necessary condition for the existence of solutions is developed which reveals that 
this dead-end point corresponds to a blow-up in the solution derivative at a point internal to the 
domain. By employing a novel asymptotic analysis in terms of a pair of small parameters, an accurate 
prediction of this dead-end point is obtained. An arc-length parameterization of the solution curve can 
be employed to continue solutions beyond the dead-end point, however, all extra solutions are found 
to be multivalued. This analysis therefore suggests the dead-end is a bifurcation point associated with 
the onset of multivalued solutions for the system. 
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1. Introduction. A micro-electromechanical systems (MEMS) capacitor consists 
of two surfaces held opposite of one another. The lower surface is a rigid inelastic ground 
plate while the upper surface is a thin elastic membrane held fixed along its boundary 
and free to deflect in the presence of a potential difference V (c.f. Fig. Il.lj) . When V is 
small enough, a stable equilibrium deflection is attained by the deflecting membrane; 
however, if V exceeds a critical value V*, called the pull-in voltage, an equilibrium 
deflection is no longer attainable and the upper surface will touchdown on the lower. 
This loss of a stable equilibrium is called the pull-in instability and the mathematical 
modeling of its onset has been the focus of many recent studies (c.f. [H Qj|] and the 
references therein for a thorough account). 

In a recent study of the physical approximations made in models of a MEMS 
capacitor (c.f. [2]), the following quasilinear, elliptic partial differential equation (PDE) 
for the dimensionless equilibrium deflection was derived: 

div ; VM = „ ^ „ , x e 0; u = 0, i6*. (l.l) 

^/TT^WW (! + «) 2 

Here, < e <C l is the aspect ratio of the device, A oc V 2 is a nonnegative dimensionless 
parameter quantifying the relative strengths of the clastic and electrostatic forces in 
the system and f2 is a bounded region in K™. The physically relevant dimensions — 
n = l, 2 — are the focus of the present work. 
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Elastic membrane at potential V 



Supported boundary 




Rigid ground plate 
Fig. 1.1. Schematic diagram of a MEMS capacitor. 



In typical applications, the aspect ratio e is a small quantity and so many MEMS 
researchers simplify (jl.ll) by linearizing the scaled mean curvature operator on the 
left-hand side, thus yielding the equation 



Au 



A 



(1 + u) 2 



x e 0; 



u = 0, x <G <9fl 



(1.2) 



This reduced equation has been extensively studied and many of its properties are well 
known (c.f. [4j [6j [19] and the references therein). One of the canonical properties is 
the existence of a critical value, A*, such that for each A < A*, (|1.2[) admits a unique 
stable solution. At the end of this branch of stable solutions there is a saddle node 
bifurcation, and accordingly no solutions of (|1.2[) exist for A > A*. In the case n = 1, 
(|1.2p has exactly two solutions for each A < A* and a unique regular solution at A = A*. 
In the case n = 2, the unstable solution branch undergoes infinitely many additional 
saddle node bifurcations which leads to higher multiplicity in the solution set (see 
the bifurcation diagrams given in Figure ITT21 where is taken to be [—1, 1] and the 



two-dimensional unit disk for (a) and|(b)| respectively). 



However, when e ^ the solution set (A, u) of (jl.ip can be markedly different from 
that of (|1.2j) . The particular focus of this paper is to catalog and analyze some of the 
profound differences between the solution structure of (jl.lj) and ()1.2|) in the singular 
limit 1 1 it | |oo — > 1- Let us begin by remarking on some differences in the n = 1 and n = 2 
cases observed from previous studies (c.f. [HED- 

When n = 1, it was shown in [I] that there exists a critical e* sa 0.35 for which the 
bifurcation curve of (jl.lj) deviates from the qualitative shape given in Fig. 1.2(a) In 



particular, if e > £*, then there exists two values A*(e) and A**(e) such that whenever 
A € (A*, A**), the stable minimal solution branch is the only solution of (jl.2j) (see 
Figure [FUjl . In the case where 17 C M 2 is the unit disk and u is radially symmetric, 
(jl.ip reduces to the ordinary differential equation 



1 



A 

(1 + u) 2 ' 



< r < 1; u'(0) = u(l) = 0, 



(1.3) 



where u(r) G (—1,0] for < r < 1. Numerical simulations have indicated the presence 
of a similar "disappearance of solutions" phenomenon, more specifically the bifurcation 
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(a)n=l, tl = {i6R; \x\ < 1} 




(b) n = 2, !i = {iel 2 : \x\ < 1} 



Fig. 1.2. Bifurcation diagrams of {1.2) . In \ (a)\ exactly two solutions are present for each 

A < A* . In \(b)\ the bifurcation curve undergoes an infinite number of fold point and hence has a 
infinite number of solutions for ranges of A. 




(a) e < e* 



(b) £ > £* 



Fig. 1.3. Bifurcation diagrams of Jl.lt for n = 1 and = {x £ R : \x\ < 1} and various e: \(a) 
for e = l/2;\m for e = 10/3. 



curve — (A(e), |u(0)|), where |w(0)| G [0, 1) — undergoes a finite number of folds before 



terminating at a single dead-end point, denoted (A*(e), |a*(e)|) (see Figures 1.4(a) 



Kb)]). Numerics show that as the bifurcation curve approaches (A*, |a*|), the derivative 
of the solution becomes unbounded at some internal point, suggesting that multivalued 
solutions may be continued beyond the dead-end point. 

Therefore to further study (|1.3p . a parametrization (r,u(r)) = (r(s),z(s)) was 
employed (c.f. [5]) in terms of an arc length s along the solution curve. The functions 
r(s), z(s) satisfy the system of ordinary differential equations: 



e 2 {z'f 



Xr' 



{l + zf 



{l + zf 



< s < L 



(1.4) 



r(0)=0, r'(0) = 1, r(£) = l, z'(0) = 0, z(£)=0, 



where r > and z > — 1 for < s < £. Note that, by the implicit function theorem, 
solutions (r(s),z(s)) of (|1.4|) can be written as z = u(r), which are solutions of (|1.3p . 
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if and only if r'(s) 7^ for all s G (0,£). Numerical simulation of (|1.4[) demonstrates 
that this arc-length parameterization allows for the bifurcation curve to be continued 
beyond (A*, |a*|) and that an infinite fold points structure, similar to that of (|1.2p for 
n = 2, is recovered (see Figures 1.4(c)| - [(d)] ). However, all of the additional solutions of 
(|1.4p beyond (A*, |a»|) are observed to be multivalued (c.f. [5]). 





(a) 



(b) Zoomed 





(c) 



(d) Zoomed 



Fig. 1.4. \(a)\ Bifurcation curves of Q1.30 computed via numerics for various e. From right 
to left the curves correspond to e = 0.05,0.1,0.5,1,2. Note that at this scale e = 0.05 and e = 0.1 
appear equal. \(b)\ A magnified portion of \(a)\ Here, the curves for e = 0.05,0.1,0.5 are seen. \(c)\ 
Bifurcation curves of 11.41 computed via numerics for various e. The dashed line represents multi- 
valued solutions of the system ljl.4[l . which naturally continue on from th e fin al classical solutions of 
11.31 1 . From right to left the curves correspond to e = 0.05,0.1,0.5,1,2. | (d)\ A magnified portion of 
\(c)\ where e = 0.05, 0.1 correspond to the right and left curves, respectively. 

The disappearance of solutions behavior is not isolated to and has arisen 

in other mean curvature type equation, notably [5], [13] and [15], which studied the 
shape of a pendant drop, electrostatic deflections of a catenoid and a one-dimensional 
Gclfand-Bratu type problem, respectively. In addition, issues of existence, uniqueness 
and multiplicity of solutions to problems of general type 

div — — = Xf(u.x), x € f2; u = 0, x £ dil 

have been a topic of recent consideration by several authors (see [HI [TT1 [HI [3] HH [12] 
[TU] 18] and the references there- in). 
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The main goal of this paper is to analyze the radially symmetric upper solution 
branch of in the limit ||it||oo — > 1 for the one- and two-dimensional unit ball. 

For the case n = 1, we analyze the upper solution branch of (jl.lj) for O = {x € 
R : \x\ < 1} using the method of matched asymptotic expansions in the limit S := 
1 + u(0) — > + . In doing so, our analysis relies heavily on the use of logarithmic 
switchback terms (cf. [HI 12). Our main result for the case n = 1 recovers the limiting 
form of the bifurcation diagram as S := 1 + u(0) — > + and is encapsulated in Principal 
Result HOI We remark that Principal Result HOI is established for any e > 0, and is 
therefore in agreement with the result (cf. [T]) that (jl.ll) admits solutions for u(0) 
arbitrarily close to —1 in the case f2 = {xeM:|x|<l}. 

In radially symmetric solutions of on the unit disc are analyzed with 

particular focus on the nature of the dead-end bifurcation. In £13 . 1 1 a rigorous necessary 
condition is established on solutions of (|1.3[) . namely that for any A < A* and e > 0, 
there exists an a*(e, A) <G R such that \u(Q)\ < \a*\ < 1. This result is proved in the 
Theorems leading up to Corollary 13 . 4.11 and demonstrates that unlike ([l.lj) in the n = 1 
case and (|1.2[) for n. = 1, 2, (|1.3[) has no solutions for u(0) arbitrarily close to —1. This 
loss of a classical solution is shown to be due to the formation of a singularity in the 
derivative at a point internal to the domain. 

To complement the aforementioned qualitative result, we employ a novel formal 
asymptotic analysis to gain insight into the disappearance of solutions at the dead-end 
point and establish a very accurate prediction of its location. The analysis demon- 
strates that the disappearance of solutions is connected in an intricate way to small 
values of the parameters e and 8 := 1 + u(0). Therefore in this case, the perturba- 
tion analysis involves two small parameters and must be performed in the distinguished 
limit e 2 /S = 5o for So =0(1). Consequently, this formal approach allows for an explicit 
characterization of the upper solution branch in terms of two functions determined by 
the solution of two associated initial value problems. These two different methods of 
computation provide a very accurate prediction of the dead-end point, as shown in 
Principal Result 13.51 by specifying a critical Sq beyond which the asymptotic solution 
fails. In agreement with Corollarv l3.4.1[ this implies that asymptotic solutions are not 
valid as S — > 0. 

In £13.31 the arc length parameterization system (|1.4j) is studied and is found to be 
amenable to the aforementioned asymptotic analysis. An accurate representation of 
the solution branch of (jl.4l) in the limit S := 1 + z(0) — > + is accordingly obtained 
in Principal Result [3.101 which is found to be valid for 5 — > 0. These additional multi- 
valued solutions therefore provide a natural continuation beyond the dead-end point 
associated with (|1.3|) . Finally, in section 21 a few open problems are discussed. 

2. One-dimensional analysis. In this section we construct an asymptotic ap- 
proximation to the maximal solution of (|1.1[) in the case where f2 is the interval [—1, 1]. 
As (jl.l[) is symmetric about the origin, we may reduce the domain to [0, 1], impose the 
condition «'(0) = and solve the resulting problem 
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in the limit 8 := 1 — |u(0)| — > + . From the bifurcation diagrams given in Figure [T3I 
we see that A = e>(l) as S —> + and, in accordance, assume 

A = V! (S)X 1 + v\ (8) Ai + . . . , (2.2a) 

where v\{8) is a gauge function to be determined satisfying v\{8) — > + as 8 — >• + . 
Furthermore, the following general form for the initial expansion of u 

u = u + vi(5)ui +vl{8)u 2 + . . . , (2.2b) 

is employed. As we will see, the final asymptotic structure of the limiting solutions as 
8 — > + is not completely obvious, and indeed our initial expansion (|2.2[) will require 
tailoring as the singularity structure of the solution manifests itself. Plugging (|2.2[) 
into (|2.ip and collecting terms, we obtain 



0, 0<x<l, ito(0) = it (l) = 0, (2.3) 

Ai 



0<x<l, «i(0) = ui(l) = 0, (2.4) 



and 



(l+e 2 K) 2 ) 3 /V (l + u ) 2 

_ A 2 (l + u )-2Aiui / 3eX«) 2 



(i + e 2 K) 2 ) 3 / 2 ; (i + w ) 3 V(i + £ 2 K) 2 ) 5/5 

u' 2 (0) = «a(l) = 0. 



, < x < 1, 



(2.5) 

In carrying out the differentiation in (|2.3p . we find that u'q = for x € (0,1), and 
therefore, uo is a linear function of x. Therefore the solution of (|2.3p that satisfies 
Uo(l) = and uo(0) = —1 is 

uq(x) = -1 + x, xg(0,l), (2.6) 

which cannot satisfy u' (0) = 0. Note that the leading order outer solution on [—1, 1] 
is uq = — 1 + \x\, which is not differentiable at the origin. This indicates the presence 
of a boundary layer in the vicinity of x = 0. 

Then using (|2.6[) in (|2.4[) - (|2.5t and solving theses equations iteratively, with respect 
to the right boundary condition, gives that 

ui(x) - -Ai(I+e 2 ) 3/2 log.T + B 1 (I-a;) (2.7) 

and 

Mx) = (1 + £2)3A? lo^ + 3A; ( l +e ^-2^A l( l +e ^ 

+ (BiA 1 (e 2 -2)Vl + e 2 - (1 + e 2 f /2 \ 2 ) log a; + d (2.8) 
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where B\ and C\ are constants pertaining to an arbitrary solution of the homogeneous 
equation and are to be determined. Therefore, from (|2.6|) - (|2.8|) . 



u = -\ + x + v 1 {8)u\(x)-\-v\{8)u< 1 + ... as(5^0 + . (2.9) 

By introducing the inner variable y = x/7, we find from (|2.9[) that the outer solution 
has the behavior u = — 1 + jy + ■ ■ ■ as x — > + . Since u = —1 + 0(5) in the inner 
layer, we set 7 = 8 and define the following inner variables: 

y = x/S, u(x) = -l + 8w(y). (2.10) 

Then substituting (|2.10[) . along with (|2.2a|) ■ into (|2.1I) we obtain 

( w' V + 



sjl + e 2 {w') 



(2.11) 



and a dominant balance requires that v\ = 8. From this we have that x = 8y, which 
implies that the local behavior of the outer expansion, (|2.9[) , is given by 



1 + SlogS (-Ml + e 2 ) 3 / 2 + (1 + ^ )3A ' ) 
8[y-X 1 (l + e 2 f' 2 logy + B 1 + (l + e 2 ) 3 A 2 



y 

3A 2 (l+e 2 ) 2 -2B 1 A 1 (l + e 2 ) 3 / 2 



2y 

+ (^JogiOCBiAiCe 2 - 2) V / l + £ 2 -(l+e 2 ) 3/2 A 2 ) 



+ 8 2 (-B lV + (Bi\i(e 2 - 2)y/l + e 2 - (1 + e 2 ) 3/2 A 2 ) logy + d) + o(8 2 ) 

as 8 — > + . However, the O(<51og(5) and C(<5 2 log5) terms cannot be matched to the 
inner expansion; so, we modify the outer expansions assumed in (|2.2[) for u and A to 
include two switchback terms: 

u = u Q + (8 log 6)u 1 /2 + 8u x + (S 2 log 8)u 3/2 + 8 2 u 2 + ..., (2.12a) 

and 

A = 8 Ai + (S 2 log (5)A 3/2 + <5 2 A 2 + . . . . (2.12b) 
By substituting (|2.12p into (|2.ip . we find that U1/2 satisfies 



1/2 



0, 0<x<l; u 1/2 (l)=0. 



Thus, the solution for u 1 / 2 is given by 

u 1/2 {x)=A 1/2 (x-l), (2.13) 
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where A]/ 2 is a constant chosen to eliminate the order O(SlogS) term. Similarly, from 
using (|2.12[) and (|2.13j) in (|2.ip we find that u 3 / 2 satisfies the ordinary differential 
equation, 



u 3/2 



2A 1/2 (l + e 2 ) 3 / 2 \ 1 



Vf+^(A 1/2 X 1 (e 2 ~2) + (l + e 2 )X 3/2 ) 
5 , < x < 1; 



«3/ 2 (l) = 0, 



whose solution is 



^4 _\ n _|_ £ 2 ) 3 / 2 

«3/a(*) = 1/2 ^ v / T+^(^i/ 2 (e 2 - 2)Ai + (1 + £ 2 )A 3/2 ) log* 



A 1/2 (l + e 2 ) 3 / 2 A 1 - ^s /2 (x-l). 



Upon substituting (|2~6)) . (|2~7|) . ([2~T3l> and (f2~T4]) into (f2~T2]) and rewriting result in 
terms of y — x/S, we have that the local behavior of the outer expansion as x — > + is 
given by 



\ V V ) 

+ s(y - Ai(l + e 2 ) 3 / 2 log y + ^ + (1 + e 2 ) 3 A?^ 

3A 2 (l + e 2 ) 2 -2 J B 1 A 1 (l+£ 2 ) 3 / 2 \ 
+ 2y J 

+ (<5 2 log 2 5) (-v/l + e 2 ^/^ 2 - 2)Ax + (1 + £ 2 )A 3/2 )) 

+ (6 2 log 6) (A 1/2 y - Vl + e 2 (A 1/2 (e 2 - 2)X 1 + (1 + e 2 )A 3/2 ) log y 

-A 1/2 (l + e 2 ) 3 / 2 X 1 + A 3/2 + B^e 2 - 2)Vl + e 2 - (1 + e 2 ) 3/2 A 2 ) 
+ S 2 (-B lV + (BiAi(e 2 - 2)^1 + e 2 - (1 + e 2 ) 3/2 A 2 ) log j/ + Ci) + o(£ 2 ). 

(2.15) 

Therefore to eliminate the O (5 log ft) terms from (|2.15j) . we choose Ax/i = — Ai(l + 
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e 2 ) 3 / 2 , and (|2.15|) becomes 

u = -l + s(y- A X (1 + e 2 ) 3/2 logy + B x + (1 + £ 2 ) 3 A 2 ^ 
V 2/ 

3A 2 (l + £ 2 ) 2 -2.B 1 A 1 (l+£ 2 ) 3 / 2 

+ 2y 

+ (<5 2 log 2 *) (A 2 (l + e 2 ) V - 2) - (1 + e 2 ) 3 / 2 A 3/2 )) 
+ (S 2 log 5) (A 1/2 y - Vl + e 2 (A 1/2 (e 2 - 2)X 1 + (1 + e 2 )A 3/2 ) logy 

- A 1/2 (l + e 2 f' 2 X 1 + A 3/2 + B.X^s 2 - 2)y/l + e 2 - (1 + £ 2 ) 3/2 A 2 
+ 5 2 (-B lV + (BiXi(e 2 - 2)y/l + s 2 - (1 + e 2 ) 3 / 2 A 2 ) logy + d) + o(<5 2 ). 

But the new order 0(8 2 log 2 J) term cannot be matched to the inner expansion, and 
again we must modify our assumed outer asymptotic expansion, (|2.12aj) . to include 
another switchback term: 

u = u Q + (8logS)u l/2 + Sui + (5 2 log 2 S)u 5/4 + (S 2 log 8)u 3/2 + 5 2 u 2 H . (2.16) 

Inserting (f2~16|) and (|2.12b|) into (f2~T|) . we obtain 
," 

l 5/4 



ii"/„=0, 0<.t<1; "5/4(1) = 0, 



which upon solving gives 

"5/4(2:) = A 5/4 (x - 1), 

where A 5 / 4 is a constant chosen to eliminate the order 0(8 2 log 2 8) term. Specifically, 
we let 



A 5/i = -^l + e 2 (A 1/2 (e 2 - 2)Aj + (1 + e 2 )A 3/2 ) 
and the local behavior of the outer solution as x — >• + is 

u = -1 + 8 (y - Ai(l + e 2 ) 3/2 logy + B 1 + (1 + e 2 ) 3 A 2 ^ 

V y 

3A 2 (l + e 2 ) 2 ~ 2giAi(l + £ 2 ) 3 / 2 
+ 2y 

+ (<5 2 log J) (-Ar(l + £ 2 ) 3 / 2 y + (A 2 (l + £ 2 ) 2 (e 2 - 2) - (1 + £ 2 ) 3 / 2 A 3/2 )) logy 

+ A 2 (l + e 2 ) 3 + A 3/2 + B^ie 2 - 2)Vl + e 2 - (1 + e 2 ) 3/2 A 2 ) 

+ S 2 (-B lV + (B1A1 (e 2 - 2)Vl + e 2 - (1 + e 2 ) 3/2 A 2 ) log y + d) + o(£ 2 ) 

(2.17) 

Besides the necessity for matching, the behavior given in (|2.17[) , along with f|2 . 10|) , 
suggests that the inner solution should be expanded as 

w = wq + (8 log S)wi/ 2 + 8wi + . . . . 
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From (|2.ip. (|2.1ip . (|2.17[) and the limiting behavior u(0) = -1 + (5 as S ->• 0+, we obtain 
the following series of inner problems: 

/ »„/ \ ' \ 1 

= -4, < y < co; w' (0) = 0, w (0) = 1, 



Vv/l + ^K) 2 )/ w o (2.18) 
io = y- Ai(l + e 2 ) 3/2 logy + Bi + o(l) as y ->• oo, 



= < y < oo; w' 1/2 (0) = 0, io 1/a (0) = 0, 

w 1/2 = -Ar(l + e 2 ) 3 / 2 y + (A 2 (l + £ 2 ) 2 ( £ 2 - 2) - (1 + £ 2 ) 3 / 2 A 3/2 )) logy 
+ X + o(l) as y -> oo, 



(2.19) 



and 



= _£, 0<y<oo; u/^O) = 0, iui(0) = 0, 



«;i = + (Bi\i(e 2 - 2)v/l + e 2 - (1 + e 2 ) 3/2 A 2 ) logy 

+ Ci + o(l) as y — >• oo, 



(2.20) 



where 



_ / y' y 2Ai 

^ : ^(l + e 2 K) 2 ) 3 / 2 J + w^- 

In particular, the solution of these inner problems uniquely determine Ai, A 3 / 2 , A 2 , -Bi, 
A 3 / 2 and C\. Note that the far field condition Wq(oo) = 1 in (|2.18[) fixes the value 
of Ai, which in turn allows the solution of (|2 . 18[) — and accordingly the value of B\ 
— to be uniquely determined. In (|2.19[) . the now fixed far field condition w' 1 ^ 2 (oo) = 

— Ai(l + e 2 ) 3 / 2 , uniquely determines the solution ^3/2,11)1/2) and, consequently, the 
value of x- By comparing (|2.19[) with (|2.17|) . the following linear equation relating the 
unknowns A 3 / 2 , x an d A 2 is obtained 



X = A 2 (l + e 2 ) 3 + A 3/2 + B 1 X 1 (e 2 - 2) v / T+7 2 - (1 + e 2 ) 3/2 A 2 . 

This process can be continued to fix the values of A 2 and C\ in f|2 .20[) . 

To determine Ai, we first multiply the ODE of (|2.18|) by w' , and observe that wq 
satisfies the following first integral for y G (0,oo): 

1 = + Ai— = —s + Ai. (2.21) 

e*y/l + e 2 {w' Q (y)) 2 Mv) £2 

Then after taking y — > 00 and using the limiting behavior of wq given in (|2.18|) , (|2.21[) 
yields 



Ai = -. 

eVl + e 2 
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To determine B\, we first solve for w' in (|2.21[) to find 



dw _ VATy/(2 - gAQw| - 2(1 - e 2 X 1 )w Q - e 2 XT 
dy e 2 Ai + (1 — e 2 Ai)uio 

An integration of this ODE yields 



< y < oo; u> (0) = 1. 



1 



e 2 X 1 + {l-e 2 X 1 )z 



Ai Jx yf{2 - e 2 Ai)z 2 - 2(1 - e 2 Ai)z - £ 2 A] 



dz = y, 



(2.22) 



where now the integral on the left-hand side can be explicitly computed and then 
expanded for y ^> 1 (see . This computation yields 



1 (!-£%) 



1 



v / ATv / 2-£ 2 Ai W ° ' VAT(2-e 2 A 1 ) 



2X^3/2 



log W 



+ 



log (A-2e 2 X 1 ) -(l-e 2 X 1 f 



O w n 



(2.23) 



VAT(2-e 2 A 1 ) 3/2 

as Wo — > oo, so that 

w ^y-X l {l+e 2 f/ 2 \ogy-X 1 {l+e 2 f' 2 {log (4 - 2£ 2 Ai) - — ^ 

as y —¥ oo. Then in comparing f|2 . 23[) with the far field behavior in (|2 . 18|) . we find 

Bi = -Ai(l + £ 2 ) 3/2 (log (4 - 2£ 2 A0 - ^ 

Next we determine A 3 / 2 - From Green's second identity and a differentiation of 
(|2~Tg)) , we obtain 



lim 

R— s-oo 



[w Jz:wi/2 — Wi/ 2 Jzw Q ) ax = hm 



fl^oo (1 + £ 2 K) 2 ) 3 / 2 



and ^fiiJo = in (0,oo), respectively. Therefore, plugging the latter into the former - 
along with the boundary conditions and limiting behavior of wq and wi — gives 



A3/2 hm / (^§) dx=-Ai, 
R-+°°Jo \w 2 J 



and consequently, A 3 / 2 = — Ai. Similarly, to determine A2, we have from Green's second 
identity, (pTO^)) and ([230|) that 



R— >oo 

Hence, 



hm / [w Q Jtwi — wiJtw Q ) dx = hm 



'(1 + £2)3/2 



R^oo (l+£ 2 K) 2 )3/2 



= A x log(4-2£ 2 A 1 ) - 



y=R 



-Si 



( 1+£ 2)3/2- 
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The determination of C\ from (|2.20[) , which is tedious and gives no special insight into 
the expansion for A, is omitted. 

We now summarize the preceding analysis for the maximal solution branch of (|2.ip . 

Principle Result 2.1. For 5 = 1 + u(0) -> + , the three-term asymptotic 
expansion for the maximal solution branch of (|2.1[) is given by 

A = SXi - (S 2 log S)X 1 + <5 2 Ai (log (4 - 2e 2 Ai) - — L= J + o(<5 2 ) (2.24) 

w/iere 

_ \/T+?-i 
e 2 v / TT£ 2 " ' 

Figure |2~T1 compares the asymptotic result of (|2.24j) and previous numerical results 
for bifurcation diagram A versus |u(0)| of (|2.1j) . As seen, the three-term result is quite 
accurate. We also remark that this asymptotic formulation does not predict the dead- 
end point A* seen in Fig. 11.31 




(a) e = 1/2 (b) e = 10/3 



Fig. 2.1. The two (dotted) and three (dashed) term asymptotic expansions from l(2~24) l of the 
maximal solution branch of 12.11 compared with the full numerics (solid). 



3. Two-dimensional analysis. 

3.1. A necessary condition for existence of solutions. In this section, we 
investigate the disappearance of classical solutions of (|1.3[) observed in [2] using tech- 
niques similar to ones introduced in [5] to study the behavior of pendant drops. For 
convenience, we introduce the change of variable p = r/e, which from (|1.3|) yields the 
nonlinear ordinary differential equation 

- ( - pU ' =|= 7 g2A v , , 0<p<£- 1 ; u'(0) = 0, uie' 1 ) = 0, (3.1a) 
with 



- 1 < u(p) < 0, < p < e 1 . 



(3.1b) 
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Here, ' now represents differentiation with respect to p. In this form, the ordinary 
differential equation in (|3.1al) admits the very advantageous geometrical interpretation: 



(p sin ip)' = 



e 2 Xp 



where ip is the angle of inclination of the solution curve u, measured counterclockwise 
from the positive p-axis to its tangent. It is important to note that these ODEs are 
equivalent on any interval in which |w/(p)| < oo. 

Now to study the non-existence of solutions of (|3.1[) . we look at a corresponding 
initial value problem, 



(p sin ip)' = 



£ 2 Ap 

(l + u) 2 ' 



p > 0; u'(0) =0, u(0) = a, 



(3.2) 



where a £ (—1,0). We begin by proving the following lemma about solutions of (13.21) . 

Lemma 3.1. If u is a solution to (|3.2p for p £ [0, a), then sin^ > for p £ (0, a), 
which implies that u is increasing on that interval. Furthermore, we have the following 
bound: 



e 2 A 



2(1 + u( P )y 



sin ib 
< < 



e 2 \ 



p 2(1 + a) 2 ' 



(3.3) 



Proof. An integration of the differential equation in (|3.2[) yields 



sin-0 



2 A 



p Jo (i + u(0Y 



d£ > 0, p£ (0, a) 



(3.4) 



Now since u is increasing on (0, a), we have a < u(p) < u(a) for p £ (0,a), and (|3.4j) 
gives ([331) . □ 

Next we prove a crucial lemma about the solutions the solutions of (|3.3p ; however, 
first we state a comparison principle [5] , which we use in our proof. 

Lemma 3.2 (Comparison Principle). Let &(£,r]) be a function such that for all 

£>o, 

Furthermore, assume that u>i(£) and w 2(C) are functions defined for £ £ [a,b] such that 



d_ 



> 



d_ 



dw 2 



on [a, b] . Suppose that 



dw x . > dw 2 



The 



Wi(b) - w 2 {b) > Wi(a) - w 2 (a). 



(3.5) 
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Moreover, equality of (|3.5[) holds if and only if Wi(£) = W2{£,) + constant on [a,b]. 
Now we can prove the following theorem. 

Lemma 3.3. Assume thatu(p;X) is a solution to (|3.2[) on an initial interval, where 
A > 0. Also, let e > be fixed. If 

(3-2-y/2) „ 4 + 3V2 
a < -1 + ^ '-£ 2 \, when < A < (3.6a) 



or 



eVA(e 2 A-8)-e 2 A-4 (3-2^2) . 4 + 3^2 . „, 

4(g 2 A + 1) <-! + " 2 £A ' wtaA> ^' ( 3 ' 6b ) 



t/ien i/iere exists a pi such that u{p) cannot be continued beyond as a solution of (|3.2p . 
Moreover, at this point (pi,u{), where U\ := u(pi), the slope of u' becomes vertical and 

2(1 + «) 2 ^ . 2(1 + M) 2 
2(1 + a) 2 



e 2 A 

where 



< ui < M < 0, (3.8) 



3a + 1 (1 + a) 2 (1 + a) Ve 4 A 2 - 12e 2 A(a + 1) + 4(a + l) 2 ,„ n , 
—\ 2?A ■ (3 - 9) 



Proof. First, we note that since u is a solution of 



simp 
P 



+ (sinip)' 



e 2 A 



(l + u)< 



(3.10) 



on the initial interval. Moreover, by Lemma 13. 11 u is increasing on this initial interval 
and we may use it as independent variable; thus, using 



-(sin 



-(cos?/)), 



dr du 

in f|3 . 10[) and then integrating the result with respect to u, we obtain 



e 2 A 



■{u - a) + (1 - cos ip(u)) > e 2 A 



1 



1 



1 + a 1 + u 



2(1 + a) 2 

where we have used (|3.3[) ; or, equivalently 

, , . e 2 X(u — a) e 2 X(u — a) 

1 - COStplu) > — r— r — r^-. 

^ v ; (l + a)(l + ti) 2(1 + a) 2 
Hence if (|3.6a|) or (|3.6bl) is true, then a vertical slope appears at a value 

Ml < M, 



(3-11) 
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where M is defined in (|3.9[) . Furthermore, the requirements in p. 21) imply that M < 0. 

First from (|3.3|) we have that u can be continued at least until p = £ 2 A/(2(l + a) 2 ), 
which implies that pi > e 2 A/(2(l + a) 2 ). 

Next, let w be defined as 

w(p) :=a + j3- ^fW~V\ P := 2 (1 + M) , 
so w(0) = a. Now, 

Wo W ')Y = £ ' 2x p < £ " x p < £2x p - tm P u >))> 

( [P ' }) (l + M) 2 (l + ^i) 2 (l + «) 2 ^ {P ' U)) ' 



for /J G (0, R), where R = minjpi, j3} and 



which by Lemma 13.21 implies 

u\p) < w'{p), u(p) < w(p), p G (0, R) 

Therefore, if /3 < pi, then w'(p) — ¥ oo as /? — > which is a contradiction, and hence, 
Pi <P. 

To get the last bound, we recall that for the solution graph (p, u) of t|3. 2|) it can be 
used as the independent variable for the solution graph (p,u). Therefore using (|3.3[) 
and fgXTJ in (|XTU)l gives 

(cos ^ > -27i £ fV (3 ' 13) 

Now let (p(u), u) be a comparison surface for u G [a, a + 2(1 + a) 2 / (e 2 A)], defined as 
(p(u),u), p(u): 



'4(1 + a) 4 



£ 4 A 2 



2(1 + a) 2 

a -' 



e 2 A 



-u , (3.14) 



with corresponding angle of inclination yj, measured counterclockwise from the positive 
p-axis to its tangent. Therefore, 

(cos ^ = ~2(iT^ 

and from ()3 . 13[) . we have (cosVOu > (cos(/j)« for each u G (a, t/), where [/ := 
min{ui,a + 2(1 + a) 2 /(e 2 A)}. Therefore, cos-0 > cosy? for it in the same interval. 
Hence the solution graph (p, u) of (|3.2[) can be continued vertically until the compari- 
son surface p,14[) becomes vertical, i.e., 

;m + af\ , 2(1 + a) 2 



e 2 A / £ 2 A 
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□ 

Therefore, we have that for a satisfying (|3.6[) . if u is a solution to (|3.2[) then its 
derivative blows- up in finite time. Furthermore, from the requirements on a in (|3.6p . 
the blow-up point {p\,u{) must happens for u\ < 0. In using this crucial fact, we can 
establish the following theorem, which rigorously proves, for all e > 0, the disappearing 
solution behavior of (|1.3p observed in [2]. 

Theorem 3.4. Let e > be fixed. Moreover, let u(r; A) be a solution to (|1.3[) , 
where u(0) = a. 

(a) If A < (4 + 3\/2)/e 2 , then 

(3-2V2) , 

a>-l + ± ^—^-e 2 A>-1. (3.15) 

(b) IfX> (4 + 3 v / 2)/e 2 , then 

eJX(e 2 X-8)-e 2 X-4 

4(^+1) ^ 

Proof. We first prove part (a). For contradiction assume that 

a < -1 + ^ ^—^-e 2 A. 

Now, since u(r) is a solution of (|1.3p . then w(p), where p = r/e, is a solution of (|3. II) . 
which in turn is a solution of (|3.2[) on an initial interval; then by the previous lemma, 
we have that u can only be continued to (pi, u(pi)), where it(pi) < 0. Therefore, since 
u is increasing on (0, pi), we obtain u(p) < u(pi) < 0, which is a contradiction, because 
u 7^ at p = e . Therefore, our assumption must be wrong, which implies that (|3.15p 
is true. 

The proof of part (b) follows similarly, except (|3.16l) is negated instead of (|3.15p . 

□ 

An immediate corollary to this theorem is the following. 

Corollary 3.4.1. Let e > be fixed and u(r;X) be a solution of (|1.3[) . Then 
there exists an a*(e,A) > —1 such that if u is a solution of (|1.3p . then u(0) > a*. 
An illustration of Theorem 13.41 is shown in Figure 13.11 

3.2. Asymptotic analysis. In this section, we use a similar analysis as in section 
[2] to analyze the upper solution branch of (|1.2p in the two-dimensional unit ball for 
e« 1. To this end, we consider the equation 

~( h TU \ J / = 7TX-v2' °<r<i; «'(o) = «(i) = o, C 

r Vvl + e 2 u' 2 J (1 + u)- 

in the limits e — > + and i^0 + , where S := 1 + u(0). In particular, the analysis will 
reveal that these two small parameters must be related together in order to facilitate 
the matching. In these limits, (|1.3p is a singular perturbation problem with an inner 
layer at r = 0; therefore, in the outer region away from r = 0, we expand u and A as 

u = u + e 2 u 1 +O{e i ), A = A + e 2 Ai + C(e 4 ), (3.17) 
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and gather terms of similar order to find 

Au ° = n A ° xa ' u "(°) = - 1 ' u o(l) = 0, (3.18a) 

(i + u Q y 

2A Ai 3A < u'$ . . 

Aui + 77T — ^3 Ul = 77T — « + on± ^2 > ui(l) = 0, 3.18b 

where A := d rr + r d r denotes the two-dimensional radial Laplacian. The general 
solution of (|3.18l) is 

u = -l + r 2 '\ A = ^; ui = ^r 2 / 3 + A sinolog r + 0) (3.19) 

for constants A, 4> — which will be determined by matching — and u> := (2\/2)/3. Note 
that Mq is not finite at r = 0, so the condition u'(0) = will need to be enforced in a 
boundary layer centered around r = 0. The value of Ai will eventually be fixed by the 
boundary condition U\(l) = 0. 

Next we analyze the boundary layer near r = by introducing the inner variables 

P = r /li w = — 1 + Sw(p) 

where 7 <C 1 is the scale of the boundary layer. Substituting these equations into (|1.3[) 
gives the following equation for w(p): 

~f -7j ^' m m ) = Sr4> o< P <co; «,(o) = i, «/(o) = o. 



A dominant balance requires that 



7 = <5 3/2 , j=5 
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where So is an 0(1) constant. We next expand w as w = wq + o(l), for 8 — >• + , and 
find that the ordinary differential equation for wq(p) is 

-(^TTTTf^) =7Z> 0<p<co; w (0) = 1, «/ (0) = 0. (3.20) 

The matching condition, from (|3. 191) . provides the leading order far field behavior for 
wq: wq ~ p 2 / 3 as p — > oo. To find the next order correction, we look for perturbations 
about this leading order form; specifically, we let wq = p 2 ^ 3 + v(p) + . . . as p — > oo, 
where i> <C p 2 ^ 3 , and retain all the linear terms to obtain the ordinary differential 
equation, 

2 An 2(5n / 

^ + 7^ + ^' = °. 

which via WKB analysis has the far field behavior 

v ~ A(Sq) sin(wlog/9 + <p(6o)), as p — > oo. 

Therefore wq ~ p 2 / 3 +« + ••• as p — >• oo, which augments p. 201) to give the full 
specification of wq as 

-f -7j ,J 0<p<oo; «, (0) = 1, «/ (0) = 0, (3.21a) 

P \ ^1 + ^oK) / w o 

u>o ~ P 2/Z + M$o) sin(wlogp + 4>(S )), as p -t oo. (3.21b) 



To carry out matching, we introduce the intermediate variable r v = r/rj(s), where 
£ 3 <C rj <C 1 as e — > + , and the corresponding order 0(e 2 ) condition 



lira - uq{ Wii ) + e 2 Ul { Wn ) + 1 - — Wo (S 3 /2 r}r v /e 3 )) ) = 0. (3.22) 
to+ £ \ "0 



fixed 

From (|3T9|) and (f3~2Tj) we have 



"ifa^) = -!r-V 2/3 rl /3 + Asin(u;log?/r t; + 
oAo 



<5o V e 3 / e 2 " 5 
as e — > + , where r„ is fixed, and (|3.22p yields 



A = ^M = a, log ^! + 

do £ 3 
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Finally applying the boundary condition Ui(l) = in (|3.19[) gives 

Mdo) .3/2^-3 



A, 



-3A - 



■sinfwlog^ £ d + ^(tf )) 



and hence, 



A = A + e 2 X 1 + 0(e 4 ) = X - 63X A(S ) sin(-\/21og(5 + 4>(S )). 



At this stage, we may fix the value of e in the main equation (|1.3[) and write So = e 2 jS 
with e 2 fixed but still 0{8). This leads to the following asymptotic result regarding 
the upper solution branch of the bifurcation curve for (|1.3[) . 

Principle Result 3.5. For solutions of (|1.3p . there is a regime where both £«1 
and S -C 1, with e 2 /S = 0(1), such that the upper solution branch of the bifurcation 
curve has the asymptotic parameterization 



| U (0)| = l-5, 



A = - — 5- A — sin 



-V21og(5 + ( — 



+ 0(S 2 ). (3.23a) 



where A(5q) and 4>{So) are functions determined by the initial value problem 



pw 



p vyi + <5oK) 2 y 

wo = P 2/3 +M s o) sin 



>o\ 



< p < oo; 



'2y/2 



\ogp- 



w {0) = 1, w' (0) = (3.23b) 
o(l) p^oo. (3.23c) 



The asymptotic parameterization (|3.23al) of the upper solution branch of (|1.3|) ap- 
pears outwardly to be defined for |w(0)| arbitrarily close to 1, potentially contradicting 
Corollary 13.4.11 However, the parameterization assumes that the quantities A{e 2 /S) 
and 4>{e 2 /5) are well defined as 5 — > + and so one can expect that A(S ) and <f>(S ) 
will not be defined for S sufficiently large. Therefore before observing the predictive 
accuracy of p.23[) . let us first consider the existence of solutions to p.23b[) . for <5 
sufficiently large. 

To do so, we follow a similar procedure outline in previous section and introduce 
a change of variables — specifically, £ = p/y/So, with wo{p) = — so that ()3.23b[) 
becomes 

-A 7= m ) =^' 2 , 0<C< TO ; «(0) = 1, v'(0)=0. (3.24) 

Hence, the mean curvature operator is isolated on the left-hand side, and (|3.24[) yields 
the geometric representation 

(esin ^' = ^v- 2 (3.25a) 
$ 9 y J 
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where ip is the angle of inclination of the solution curve. Noting that (|3.24[) and (|3.25a[) 
are equivalent on any interval in which is bounded, we look at (|3.25a[) . coupled 
with the initial condition 

v(0) = 1, (3.25b) 

to study the nonexistence of solutions of (|3.24[) . Also note that if v satisfies (|3.25l) . 
then the condition v'(0) — is redundant, which can be seen by integrating (|3.25a[) 
and then taking £ — >• + . 

Lemma 3.6. If v satisfies (|3.25j) on [0,a), then simp > for £ £ (0,a) ; which 
implies that v is increasing on that interval. Furthermore, we have the following bound: 

2S < sin|(0 < 2|o_ (3 ^ 6) 



9v(o 2 e 9 

Proof. Integrating (|3.25ap yields 

sin^o = ^-jf v «fa)~ 2 d ^>°- 

for all £ G (0, a), and the results follow as in Lemma I3TT1 □ 

Now, we can prove the main lemma which leads to our desired main result for the 
nonexistence of solutions of (|3.23b|) . 

Lemma 3.7. Assume that v is a solution of (|3.25a[) on an initial interval with 
v(0) = 1. If 

- 9(2^+3) 
Oo > o ■= 2 ' {3.27} 

then there exists a value £i in which cannot be continued beyond as a solution of 
(|3.25[) . Furthermore, at this point where V\ '■= u(£i), the slope of the solution 

curve is vertical and the following bounds hold: 



9 9 66 - 9 -y/46*- 1086 + 8l \ 
2^ <6 -2^ 4^ ' (3 ' 28) 



_9 < 6* -9- V ^ 108^o + 81 

2o 4d 



Proof. By Lemma 13.61 wc know that on an initial segment v is increasing and 
consequently, may use it v as an independent variable; thus, (|3.25p gives 

— (cosV') u = -g-u , (3.30) 

Then integrating with respect to v and using (|3.26p yields 

AS ( 1\ 2S a 
1-cos^H > — (l--j -—(«-!), 



ANALYSIS OF PRESCRIBED MEAN CURVATURE EQUATION 



21 



or 

2\ 

H — ) > cosib(v); 
vj 

hence if (|3 . 2T[) is true, then a vertical slope exists at a value V\ where 

6<y - 9 - JA51 - 108(5 + 81 

Vl < AX • 

4d 

First from (J3T26J) , we see that v can be continued to at least £ = 9/(2<5o), which 
implies that £i > 9/(2<5o). 

Next, we compare the solution to 

W(Q:=l + e-V^, g : = ^. ^-»-^- -i°g*±a)'. 
Note that 

for £ e (0,min{fi,/3}) and W(0) = 1. Therefore, by LemmaHU 

< W(0 < 1/(0. C € (0, minfe, /?})• 

Now if p < then < «'(£), for £ € (0,/3), which implies that i/(/3) = oo and 

yields a contradiction. Therefore, we must have £i < (3. 

To get the last bound we note that for the solution graph (£, v) of p.25[) v may be 
used as the independent variable. Hence using (|3.26[) in (|3.30l) gives 



(cosi/j) v > — 




(3.31) 



Then in introducing a comparison surface 



(3.32) 

for v G (1,1 + 9/ (2 Jo)) with corresponding angle of inclination </?, we obtain (cos (p) v = 
-9/(2<5 ); thus, from (l3~3T1) 

(cOSf/>)i> > (cOS(/?)„ 

for each i> <E (0,min{wi, 1 + 9/(2<5o)}), which implies cos'0 > cosy for each v in that 
interval, and therefore, the solution (£, v) of (|3.25[) can be continued vertically until 
the comparison surface (|3.32p becomes vertical. That is, 

vi > u;(9/(2<5o)) = l + 9/(2<5 ). 







(v + 


2<5o 



(Z(v),v), Z(v)-=\IOk- 



81 



2S 



22 



N. D. BRUBAKER AND A. E. LINDSAY 



Therefore, since (|3.25[) and (|3.24p are equivalent on (0,£i), then the derivative of 
(13.241) also blows-up in finite time, which after changing variables back to p yields the 
following result concerning (|3.23b[) . 

Theorem 3.8. There exists a value Sq such that for Sq > Sq no solutions of 
(I3.23bp exist. Furthermore, 



S* Q < S := 9(2 ^ +3) « 26.2279. 



Remark. An integration of the (|3.23bp gives Sq sa 18.142468. indicating that So is not 
a particularly tight upper bound on Sq. 

As a result of this theorem, an expansion for the dead-end point (A*(e), |a*(e)|) can 
now be extracted from (|3.23a|) ; specifically, since (|3.23b[) has no solutions for So > Sq, 
the asymptotic approximation (|3.23a[) fails at e 2 /5 = Sq, or S = e 2 /Sq. Therefore, using 
these values in p.23aj) . the following asymptotic result for the dead-point of (| 1 . 3[) is 
established. 



Principle Result 3.9. For e <C 1, the dead-end point of the upper solution 
branch of the bifurcation curve of (|1.3j) has the asymptotic expansion 



\a*(e)\ = l- c - + 0(e 4 ), 



A*(,)^4^sin 



9 3 S* 



-V2\o g - + 4>(6Z) 
°o 



(3.33) 

+ 0{e A ). 



We remark that the ability of Principal Result 13.91 to predict the dead-end point 
associated with radially symmetric solutions of (jl.ip in the case n = 2, is not matched 
by the asymptotic analysis leading to Principal Result 12.11 for the n = 1 case. This 
discrepancy suggests that the dead-end phenomena exhibited in both the n = 1 and 
n = 2 are not qualitatively similar. 



In order to study the quantitative accuracy of Principal Result 13.91 it is necessary 
to obtain the functions A(Sq) and 4>(Sq), which are readily acquired by solving (|3.23b[) 
numerically, then subtracting off the growth term p 2 / 3 and applying a least squares 
fit to the remainder (see Figure |3~2|) . In Figure [3~3l comparisons of the full numerical 
solution of the upper branch of the bifurcation curve and asymptotic prediction of 
(|3.23a[) are displayed; furthermore, the agreement is observed to be very good. 
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(a) A(S ) 



(b) <«<5o) 



Fig. 3.2. Graphs of A(8q) and <j>(So) against 6q. The numerical integration fails abruptly at 
roughly <5 = <5q ~ 18.142468. 





(a) £ 2 



0.025 



(b) a* 



0.1 



Fig. 3.3. Comparison of the full numerical solution of the bifurcation curve for Hl,3j l (solid) with 
asymptotic formula Il3,23al) (dashed) and the solution curve of il.3t for 8 = (dash dot). 



In Figure kOI a comparison of the numerical and asymptotic values for the location 
of the dead-end point is shown; note that agreement is very good, as in each case the 
asymptotic error is C(e 4 ). 

In Figure 3.5(a) the numerical (solid) and global asymptotic (dashed) solutions of 
(11.3[) at the dead-end point for e 2 = 0.2 are displayed. As expected (see section [3~Tj) . 



the tangent of the solution curve is almost vertical, indicating that the derivative of 
the solution is becoming unbounded. Indeed, when solutions w' (p;6o) of (|3.23b[) arc 
plotted for several So a blow-up in w' (p; So) as So — > Sq is observed (see Figure 3.5(b) |. 



This suggests that beyond the dead-end point solutions of (|1.3p cannot be represented 
by a function of a single variable. Therefore in the next section (|1.3[) is re-parameterized 
in terms of arc length along the solution curve s, and consequently becomes a system of 
coupled ODES. An asymptotic study of this coupled system reveals that multivalued 
solutions of ([l.ip are present beyond the dead-end point of the bifurcation diagram. 
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(a) e 2 /8* =1-K(e)| 



(b) A.(e) 



Fig. 3.4. Comparison of the asymptotic prediction, 13.331 1 . (dashed line) of the dead-end point 
(A*(e), |a*(e)|) with full numeric computations (solid) for: \(a)\ the 0(e 2 ) correction of \a m (e) | ; | (b)\ 
A*(e). Notice that the scale on the y-axis of the right figure is quite fine and so the agreement for 
A* (e) is in fact better than the figures makes it appear. 




(a) Global Approximation - e 2 



0.2. 



(b) Blow up of u>' (p; So) as <5q — > 6J. 



Fig. 3.5. \(a)\ The numerical (solid) and global asymptotic (dashed) profile in the boundary layer. 
The curve is almost vertical indicating that the derivative is becoming infinite. \(b)\ The derivative 
function w' (p;So) plotted for several <5q . As So — > Sq, it observed that w' (p;So) appears to develop a 
singularity at a finite p* . 



3.3. Arc length asymptotic analysis. In this section, we analyze parameter- 
ized solutions (r(s), u(r(s))) = (r(s),z(s)) of (|1.3|) where the arc length s satisfies the 
relationship dr 2 + e 2 dz 2 = ds 2 . When reformulated in terms of r(s) and z(s), (|1.3|) 
becomes a coupled set of ODEs: 



e 2 {z'f 



(l + z) 2 



z" + 



Xi 



r (l + z) 2 ' 



< s < 



r(0) = 0, r'(0) = 1, r{t) = 1, z'(0) = 0, z(£) = 0, 



where < e 2 <C 1 and A and £ are unknown parameters to be determined. Furthermore, 
to facilitate the analysis of the upper solution branch, we impose the condition z(0) = 
— 1+5 and study (|1.4|) in the limits e — > + and 6 — > + , where the relationship between 
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these two small parameters is to be determined. In the outer region away from s = 0, 
we expand r, z, A and £ as 



(3.34a) 



r(s; e) = r {s) + e 2 n{s) + e 4 r 2 (s) + 0(e e ), 
z{s; e) = z (s) + e 2 zi(s) + e 4 z 2 (s) + 0(e 6 ), 

and 

X(e) = X +e 2 \ 1 +e i X 2 +O{e 6 ), t(e) = £ + e 2 £ x + e A l 2 + G(e 6 ), (3.34b) 
which upon substituting into (|1.4[) gives 



r £ = , z» + ?M = M. 0<S <1; 

° r (l + z ) 2 (3.35) 

r (0) - 0, r' (0) = 1, r (£ ) = 1, z(0) = -1, z' Q (0) = 0, z (£ ) = 0, 

at order 0(1). Therefore in solving (|3.35l) we find 

r (s) = s, z (s) = -1 + s 2 / 3 , A = ^, e = l. (3.36) 

However, z' (0) ^ and therefore, we have a boundary layer at s = for z(s). Next, 
from (HHI), (|3"j&|) and (|3~3B)) we have 

n(o) = 0, ri(0) = 0, n(l) = 2 X (1) = 



at order C(e ). The solution for r±(s) is 



2 « A . 2 1/3 



n(s)=^--£i)s--s L '- s (3.38) 

where the condition r[(0) = will be enforced in the boundary layer at s = 0. Then 
using ()3.38j) in (|3.37[) . we deduce 

which upon solving gives 

aw = (27Al+ 3 8 6 ' 12fl) g2/3 ~t +Ai sin ( w log s + ^ ' ^ 3 - 39 ) 

where u := 2^/2/3. Here, Ai and y>i are constants that will be determined by matching 
and the value of Ai will be determined later by applying the condition Zi(l) = —24/3. 
In order to fix the value of £\, an expansion to higher order is required. Accordingly, 
we use (|3.34j) in (|1.4|) to find a system of ODEs at order C(e 4 ) (see[C]), which upon 
solving gives 

Ki 64-4 + 36^5^0! 1/3 
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where 



and 



where 



2 + 4^2 Ai cos (lu log s + 0i ) - 16 A sin (u log s + fa) 



27 

, 2 2^i 9 2 20 

A' 2 := ( 27 - — + e i - h - C 2 - -uAi cos fa - —A sin0i 



Z2 ( s ) = ~2j3 + giJ5 + \ Al sin<?il + ^2 sin (w logs + fa) + K 5 s 2/3 , 



2 A\ 58^ A . , , . 20 . 

:= 81 + ~2 81 - cos(wlogs + 0i)- —A sin (a; logs + <?!>i) 

+ ^A 2 cos(2a;logs + 20 1 ) + ^^-A\ sin(2a; logs + 200, 
38 19 



4 - 364 + 54£ 2 - 544 - 54C 2 - 8 V2A t cos 0i - 40A sin fa 

Ka . = , 

18 

484 - 36^ - 16 - 162A 2 + 243A 2 + 108C 2 

K5 '- = 324 

9 A 2 cos 20i - 2Ai (34 - 2) sin 0i 
+ 18 ■ 

Next we introduce the inner variable p = s/7, which after plugging into zq gives 
the near field behavior z = — 1 + 7 2 / 3 p + ... as s 4 + . Moreover, z = — 1 + 0{5) in 
the inner layer, and as a result, we choose 7 = <5 3 / 2 . Then for matching we write the 
outer solution, (|3.34aj) . in terms of the inner variable, p = s/<5 3 / 2 . to obtain 



S 3/2 (p-^P 1/3 + f^)+^2 



+ ^ (So (§ ~h)p + ^(^-4+^1^01)^ + 0(s7/2) 



(3.40a) 



2 = -1 + <5 ^p 2/3 - + £ Aisin (wlogp + v^logtT 
+ ^((^^^-^isin0i),V3 



(3.40b) 



2 , . \ 

-A sin 0! + (5 2 A> sin (w log (<5 3/2 p) + 2 ) + 0(S 3 ) 



3 

as s — > + . It will be seen that the inner solution cannot be matched to the order 
0(S 2 ) of (|3.40a[) . For this reason, the constant C 2 is chosen such that K2 vanishes, 
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which causes K± to vanish and gives the reduced local behavior 
2S ° -1/3 M n( n -V*\ 



J3/2 



+ S (So (I - 4) p + ^(^-4 + 36^^^)^ 



0(S 2 ) 



(3.41a) 



-1 + 5 



p 2/s _ f_o + ^ gin ^ bg p + V2 log + (^1) + 0(p- 2 / 3 ) 
+ S ( ( (4 '^ l)(5 ° - M! sin fc) p 2 /3 + (1)\+0(S 2 ) 



(3.41b) 



as s — > + . As a consequence, these local expansions motivate us to introduce the 
following local variables within the vicinity of s = 0, 



p = sl^l\ r(s)=8^ 2 R(p), z(s) = -l + SZ(p), 
which transform (11.41) into 



R" = 



\Z' t {Z'f 



R 



R'Z' _ XR' 



R 



(3.42) 



(3.43) 



Here the dominant searings require e 2 /S = 8q, where 8q = 0(1). Therefore, expanding 
R, Z and A as 

R = Ra + SR^+O^S 2 ), Z = Z a + 8Z1 + 0{8 2 ), A = A + 880M + 0{8 2 ), (3.44) 
respectively, we find that the leading order problem for the inner solution is 



oil _ c ^0 Z , c C^O 



Ro 



7 2 

A) 



< p < 00; 



(3.45) 



i?o(0) = 0, R' o (0) = 1, Z o (0) - 1, 4(0) = 0. 



To find the far field behavior of i? and Z , we assume Ro^p + V&sp— > 00, where 
F < p, and Z ~ p 2 / 3 + VK as p — >• 00, where W <C p 2 / 3 . Substituting these relations 
into (I3.45|) . gives asymptotic differential equations for V(p) and W(p), 



W 2A , 



27 



32£ 1 

"sTp 2 



as p — > 00, 



whose solution is V ~ — (2<5o/3)p 1 / 3 , ~ — <5oAo + Ai(So) sin(w log p + 4>i(So)) as 
p — > 00. Hence, the far field behavior for the solution, Rq and Zq of (|3.45p is 



Ro(p) = p-^p 1/3 + o(1), 



MP) = P 2/3 -^oA + ii(5 )sin (wlogp + ^i(* )) +o(l), 



as p — > 00. (3.46) 
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Proceeding to 0{8) terms, we substitute (|3.44[) into (|3.43|) and collect the 6 terms to 
find that R\ and Z\ satisfy 



R" = So 



7 2 
A) 



R 2 



, < p < oo, 



ryit _ ^oAi-Ro + Ao-Rj 2\ ZiR' RiR' Z' R[Zq + R' Z[ (3.47) 

1 — "rpl ~r7?> ' o2 TJ ' ^ 
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i?i(0) = 0, i?i(0) = 0, Zi(0) = 0, Z((0) = 0. 



From (|3.41[) we expect the far field behavior of both R\ and Z\ to grow algebraically. 
Therefore, we assume R\ ~ ap Q and Z\ ~ as p — > oo and substitute this behavior 
into (|3.47p along with the far field behavior of Rq and Zq. After a dominant balance, 
this yields 



aa(a — l)p° 



2£q Al 5/3 



6)808 - I)/" 2 ~ <5 A 1/9 - 4/3 - 6 (/? + 2A ) /" 2 , 
as /? — s- oo . Consequently, 

a = -d Ai, « = -, 6= — - — , 



which implies that 

Ri ~ -6l\ip 1/3 , Z x 



3^qAi 2/3 

A P ' 



as p — > oo. 



(3.48) 



As a result, (J3T44J) , (|3.46j) and (|3.48j) give the following far field behavior of the inner 
solution: 



z = 



p 2/3 - S X + Ii (5 ) sin L log p + fa (<5 )) + . . .) 

-0{5 2 ) 



(3.49) 



/ 3£ Ai 2/3 

V 4 P 



as <5 — > + and p — > oo. Then to match we compare (|3 .42[) . using (|3.49[) . with (I3.4ip to 
get 



ii(^o) 
<5o ' 



ii(*o)~ V21og5, Ai = --A lS in^. (3.50) 



Note that the boundary condition zi(l) = — 2£i/3 is automatically satisfied by the 
value of Ai determined in (|3.50[) . By returning to the definition of A made in (|3.34bl) 
and recalling that e 2 /S = 5$ , a two term expansion of A is now given by 



A = 5 - Ai 

9 3 



V21og(5 
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Next, we fix e in the governing equation, (|1.4[) ; therefore, for our asymptotic analysis 
to remain valid, we need e 2 /S = Sq = 0(1), with e fixed, which leads to the following 
asymptotic result regarding the upper solution branch of the bifurcation diagram of 

cm. 

Principle Result 3.10. For solutions of l|1.4p . there is a regime where both 
e <C 1 and i « 1, with e 2 /S = 0(1), such that the upper solution branch of the 
bifurcation curve has the asymptotic parameterization, (X(S;s), |z(0)|), where 



\z(0)\ = l-6, A = A -<^ii(% Isin 



i(-) -V21og5 



0{8 2 ). (3.51) 



Moreover, Ao = 4/9, and A\(So) and 4>i(So) are functions determined by the far field 
behavior of Zq , 



p 2/3 _ ^0 + ^ 1 ( ( J ) sin Lj\ og p + fa(S ) 



9 



Zo(p) 

of the initial value problem 



9 Z 2 Rq 



Z n , R 'o z o _ 4 

Ro 9 Z 2 ' 



i?o(0) = 0, R' o (0) = 1, Z o (0) = 1, Z' o (0) = 0, 



< p < oo; 



(3.52a) 



(3.52b) 



where 



2(5 

Ro(p) = P~ -^p 1/3 + o(l) as p^oo. 



(3.52c) 



To study the accuracy of this result, we again need to compute the functions A\ (Sq) 
and 4>i (Sq), which is done by subtracting the growth term and applying a least squares 
fit to the remainder. The graphs of these functions are displayed in Figure 13.61 As 
expected, these new functions are continuations of the old functions found in Principle 
Result [331 

A combination of the asymptotic formula (|3.51[) and the numerically obtained 
functions Ai(Sq) and </>i(<5o) allow for a reconstruction of the bifurcation diagram (see 
Figure [377| . Numerically, Ai(5q) appears to grow linearly as So — > oo, which would 
indicate that SAi(e 2 /8) is finite as S —> 0. Therefore the analysis predicts that the 
upper solution branch of (|1.4[) undergoes infinitely many fold points in a way similar 
to the upper branch of (|1.2j) for the two-dimensional unit disk. This prompts the 
following conjecture. 

Conjecture 3.11. For e > fixed and sufficiently small, the upper solution 
branch of the the bifurcation diagram of (|1.4j) undergoes infinitely many folds and as 
1 2 (0)| — > 1~ , A > goes to a finite value that is bounded away from zero. 

In Figure 13.71 (|3.5ip is compared with the numerically computed bifurcation dia- 
gram of (|1.4[) . From this we see that the observed agreement is very good. 

4. Conclusion. In this work, we have analyzed the upper branch of solutions to 
equation (jl.lj) in the limit as |u(0)| —> 1~ with particular focus on one and two spatial 
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Fig. 3.6. Graphs of Ai(<5o) and </>i(5o) against 80 computed from 113. 5211 . The solid line indicates 
where the coefficients agree with those compute from J3- 23D - 




Fig. 3.7. Comparison of the full numer ical solution o f the bifurcation curve for l|1.3|l (so lid) 
with asymptotic formula ( I3.23atl (dashed) for: \ (a)\ e = 0.05; {(b)\ s = 0.1. Note that in both \(aj\ and 
(b)\ the asymptotic prediction agrees extremely well with the full numerical solution. 



dimensions. In both cases, there are marked differences between the solution structure 
for e = and e > 0, specifically a disappearance of solutions is observed, i.e., (jl.ljl 
does not necessarily have a solution for all |u(0)| G [0, 1). 

For the case n = 1, there are two qualitatively different bifurcations structures 
associated with which are separated by the cases e < e* and e > e* (c.f. Fig. II. 3|) 
for a critical value e* ~ 0.35. In the case where e > e* , the upper solution branch of 
equation (|1.1[) is defined only for |u(0)| € [0, a**(e)) U (a*(e), 1), where < < 
a*(e) < 1, while for e < e* equation (jl.lj) has a solution for any |u(0)| G [0,1). 
Therefore, in the case n = 1, equation (jl.lj) has solutions for |u(0)| arbitrarily close to 
1 , and these solutions have been constructed using singular perturbation techniques in 
the limit |u(0)| — > 1 _ . The limiting form of the bifurcation diagram, encapsulated in 
Principal Result [270 has been found to be highly accurate, however, in contrast to the 
asymptotic analysis of symmetric solutions of (jl.ljl in n = 2, it does not predict the 
dead-end point (A*(e), a*(e)). 
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In the case n = 2, we have shown that solutions of (|1.3[) do not exist for |u(0)| 
arbitrarily close to 1 for any e > 0. It is observed that as the upper solution branch 
is traversed, a singularity in the first derivative of the solution develops interior to 
the domain and at this singularity, the branch of solutions ends abruptly at a single 
dead-end point. Our asymptotic analysis allows for an accurate prediction of this point 
to be made by relating it to a singularity in an associated initial value problem. The 
analysis predicts that the singularity occurs for a fixed value of e 2 /(l + |u(0)|) and 
therefore establishes a relationship between a given e and the dead-end point. In each 
case, the asymptotic parametcrizations obtained for the solution branch compare very 
well with full numerical solutions. The ability of the asymptotic analysis to predict the 
dead-end point when n = 2 case but not when n = 1, suggests a different underlying 
mechanism is responsible for the phenomena in each case. 

Finally, by employing an arc length parameterization of solutions to (|1.3|) , we find 
and analyze a new family of solutions emanating from the dead-end point. These solu- 
tions are found to be multi- valued and provide a natural continuation of the bifurcation 
curve beyond the dead-end point which retains the infinite fold points feature of the 
e = problem. 

The main limitation of our study is that for n = 2, we deal only with radially 
symmetric domains. Though our analysis has revealed interesting structure associated 
with solutions of (|1.3[) . an investigation of on more general domains would be 

desirable; specifically, can a result like Theorem 13.41 be formulated, where |m(0)| is 
replaced by ||u||oo- 

Additionally, it would be interesting to study solutions of in higher spatial 

dimensions. Rigorously it has been shown that when fl is the n-dimensional unit 
ball and 2 < n < 7, the bifurcation diagram of (|1.2j) exhibits the infinite fold points 
structure [7] . What then is the effect of positive e on the bifurcation structure of (jl.ip 
when n > 3? Another interesting avenue for future investigation is the dynamic version 
of (jl.ip . namely the equation 

u t = div — = - — -k, x S S2; 

v/l + e 2 |Vw| 2 (l + u) 2 (4.1) 

u(x,0) = 0, u = 0, xedQ,. 

Is there an equivalent of disappearance of solutions for (|4.ip . i.e., does u or its deriva- 
tives exhibit a singularity at some finite t before reaching u — —11 

Acknowelegments. N.D.B. would like to thank J. A. Pelesko for many useful 
discussions. 

Appendix A. Here we compute the integral given in (|2.22[) and then expand for 
y 1. First, we split the aforementioned integral into two parts: 

Woiv) e 2 \ 1 + {l-e 2 \ 1 )z T T ^ 

v ; = dz = h + 1 2 (A.l) 

y/{2 - e 2 Ai)z 2 - 2(1 - e 2 Xi)z - e 2 Ai 
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where 



h := s 2 Xi / ; dz, 

y/{2 - e 2 Ai)z 2 - 2(1 - £ 2 Ai)z - £ 2 A X 

rwo(y) 

h ■= (1 - e a Ai) / , dz. 

y/(2 - £ 2 Ax)z 2 - 2(1 - e 2 Ai)z - s 2 X 1 



In computing I\ and I2, respectively, we obtain 

v /((2-e 2 A 1 )z-(l-e 2 A 1 )) 2 -l 

£ 2 Ai /■(2-e 2 A 1 )«;o-(l-e 2 Ai) 1 



fwa(y) 1 

h ■= e 2 \i\/2-e 2 \ 1 / = dz 

h . do _ c 2i, w _ n _ c2\\\ 2 _ 1 



V2 - e 2 Ax A VF^T 
e 2 A! 



cosh -1 ((2 - e 2 Ai)u; - (1 - £ 2 Ai)) 



and 



V2 - e 2 Ax 



J a := (l-e 2 A 1 )v / 2-e 2 Ai y 



v /((2-e 2 A 1 )z-(l-e 2 A 1 )) 2 -l 

(1 - £ 2 Ax) /•(2-e 2 A 1 )^o-(l-e 2 A 1 ) £ + (1 _ ^ 



di 



(2-£ 2 A 1 )3/ 2 i 1 v ^2— 1 

v /((2-£ 2 A 1 ) Wo -(l-£ 2 A 1 )) 2 -l 

cosh" 1 ((2 - e 2 \i)w - (1 - £ 2 Ai)) 



(2-e 2 Ai) 3 / 2 

(l-£ 2 Ar) 2 



(2-£ 2 Ai) 3 / 2 

Therefore, since t/»L«)o>l, which implies 

cosh" 1 ((2-£ 2 A!) Wo - (1 -£ 2 Ai)) = logu> +log(4-2e 2 Ai) + O (w" 1 ) 
^/((2-£ 2 A 1 ) U ; -(l-£ 2 A 1 )) 2 -l = (2 - e 2 AOt«o - (1 - £ 2 Ai) + O (w^ 1 ) 



and hence, 

r £ 2 Ai . £ 2 A 1 log(4-2£ 2 A 1 ) _ x 

J 1 = — log W)n H , + O ( U> n ) , 

V2 - £ 2 Ai 8 V2 - £ 2 Ar V ; 

(l-£ 2 Ai) (1-£ 2 A0 2 , , kn . 

(l-£ 2 A 1 ) 2 log(4-2£ 2 A 1 )-(l-£ 2 A 1 ) 2 _ 

+ ( 2 - £2Al ) 3 / 2 + ° K } 

as u>o — > 00. Therefore, from (|A.1[) and (|A.2[) we obtain 

(l-£ 2 Ar) 1 , log(4-2£ 2 A 1 )-(l-£ 2 A 1 ) 2 _ 
7 = w o H ^ log w H T75 h C ( w ) • 

V2 - £ 2 Ai (2-£ 2 Ai) 3/2 (2-£ 2 Ai) 3/2 
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Appendix B. To solve (1 1 .4|) numerically, we use a shooting method. That is, we 
impose the initial conditions 



r(0) = 0, r'(0) = £, z'(0) = 0, u(0) = a, 



(B.l) 



where a G (-1, 0) and find (A, £) such that F(X, I) := [r(l; A, ^) — 1 z(l; A, J?)] T = 0. 
To do so, we apply Newton's method and iterate as 



A n +i — A„ — 



dr dr 
— (1;A„,C) — (l;A„,£ n ) 

dz dz 



r(l;X n ,£„) - 1 
(1;A„,£„) 



where A„ = [A„ £„] . Therefore at each step we need to find r\(l; A„, £ n ), rg(l; A„, £ n ), 
Z\(l; X n ,£ n ) and z e (l; X n ,£ n ). To this end, we differentiate the ode given in (| 1 .4[) and 
the initial conditions (|B.1[) with respect to A and separately with respect to £ to get 
two auxiliary problem for (r\, z\) and (r^ , zi), whose solutions evaluated at £ = 1 yield 
our desired result. 

Appendix C. Here are the order 0(e 4 ) outer ODEs for (jl.4l) : 

„ Mi M 2 

r 2 s 7/3 + s 5/3 ' 



r 2 (l)=4(4-A )-4, 
„ , 1 / , 2A 2r 2 
s s z 



2r'« 



Ms t M 4 M 5 

s 4/3 s 8/3 



4(154-8) 2^ 2Ai£ 
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— I — (sm0i — v 2 cos 0i J 



(C.l) 



where 



Mi 
M 2 
M 3 
M 4 

Mr 



144^^1 cos log s + 0i) + 144Ai sin log s + 0i) + 8 



243 

8- 124 - 72Ai sin 0i 
243 

32(3-24-9^1 sin0i) 
729 

227 + 484 - 36^ - 36Ai(3£i - 2) sin0i - 324A 2 sin 2 0i 



4 (4 + 54 v / 2Ai cos log 



243 

s + 01 



I8OA1 sin ( 2^2 log s + 0i 



4AI sin 2 logs 



729 



and we have simplified the result using (|3.36[) . (|3.38[) . (|3.39[) and (|3.50D . 
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